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The  morphologically  similar  taxa  Anopheles  calderoni,  Anopheles  punctimacula.  Anopheles  malefactor  and 
Anopheles  guarao  are  commonly  misidentified.  Isofamilies  collected  in  Valle  de  Cauca,  Colombia,  showed  morphologi¬ 
cal  characters  most  similar  to  An.  calderoni,  a  species  which  has  never  previously  been  reported  in  Colombia.  Although 
discontinuity  of  the  postsubcostal  pale  spots  on  the  costa  (C)  and  first  radial  (R)  wing  veins  is  purportedly  diagnostic 
for  An.  calderoni,  the  degree  of  overlap  of  the  distal  postsubcostal  spot  on  C  and  Rt  were  variable  in  Colombian  speci¬ 
mens  (0.003-0.024).  In  addition,  in  98.2%  of  larvae,  seta  1-X  was  located  off  the  saddle  and  seta  3-C  had  4-7  branches 
in  86. 7%  of specimens  examined.  Correlation  of DNA  sequences  of  the  second  internal  transcribed  spacer  and  mtDNA 
cytochrome  c  oxidase  subunit  I  gene  (COI)  barcodes  (658  bp  of  the  COI  gene)  generated  from  Colombian  progeny  ma¬ 
terial  and  wild-caught  mosquitoes  from  Ecuador  with  those  from  the  Peruvian  type  series  of  An.  calderoni  confirmed 
new  country  records.  DNA  barcodes  generated  for  the  closely  related  taxa,  An.  malefactor  and  An.  punctimacula  are 
also  presented  for  the  first  time.  Examination  of  museum  specimens  at  the  University  of  the  Valle,  Colombia,  revealed 
the  presence  of  An.  calderoni  in  inland  localities  across  Colombia  and  at  elevations  up  to  1113  m. 
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The  Central  and  Southern  American  mosquito  spe¬ 
cies  Anopheles  (Anopheles)  calderoni  Wilkerson,  1991, 
Anopheles  (Anopheles)  guarao  Anduze  and  Capdevielle 
(1949),  Anopheles  (Anopheles)  malefactor  Dyar  &  Knab 
1907  and  Anopheles  (Anopheles)  punctimacula  Dyar 
&  Knab  1906  (Arribalzagia  Series,  Laticorn  Section) 
(Harbach  2004)  are  highly  isomorphic  and  closely  re¬ 
lated  (Wilkerson  1990,  1991).  Following  reassessment 
of  specimens  from  both  the  type  series  and  type  local¬ 
ity  [Panama,  Colon  (Canal  Zone)]  of  An.  punctimacula, 
Wilkerson  (1990)  detected  exemplars  of  An.  malefac¬ 
tor  in  the  syntype  series  and  resurrected  the  nominal 
species  from  synonymy  with  the  former.  According  to 
Wilkerson  (1990,  1991),  An.  punctimacula  is  reported 
from  Mexico  to  northern  Colombia  and  Venezuela, 
whereas  An.  malefactor  appears  to  be  limited  to  Pana¬ 
ma  and  northeastern  Colombia.  So  far,  the  distribution 
of  An.  guarao  is  thought  to  be  limited  to  its  type  local¬ 
ity  of  Caripito,  Estado  Monagas,  Venezuela  (Anduze  & 
Capdevielle  1949). 
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Further  assessment  of  An.  punctimacula  s.l.  specimens 
from  Piura,  Salitral,  Peru,  revealed  yet  another  cryptic 
species,  which  Wilkerson  (1991)  formally  named  An. 
calderoni.  The  description  suggested  that  prior  records 
of  An.  punctimacula  in  the  low-lying  areas  (<  250  m)  of 
the  Pacific  Coast  of  Colombia,  Ecuador  and  Peru  (Levi- 
Castillo  1949,  Russell  et  al.  1963,  Calderon  et  al.  1974) 
could,  at  least  in  part,  refer  to  An.  calderoni.  Based  on 
characters  of  the  adult  females,  Rubio-Palis  and  Moreno 
(2003)  confirmed  An.  calderoni  from  the  states  of  Tachira 
and  Barinas,  in  Venezuela.  In  light  of  this  geographi¬ 
cal  discordance,  its  presence  in  Colombia  and  Ecuador 
seems  highly  likely.  That  An.  calderoni  has  been  found 
positive  for  Plasmodium  vivax  in  the  Peruvian  regions  of 
Lima,  Salitral,  Querecotillo,  Lourdes  Paz  and  Palo  Blanco 
(Kroeger  &  Alarcon  1993,  Calderon  et  al.  1995),  makes 
this  taxonomic  query  also  one  of  biomedical  importance. 

Literature  records  and  museum  holdings  indicate 
only  the  presence  of  An.  punctimacula  and  An.  male¬ 
factor  in  Colombia  (Gonzalez  &  Carrejo  2007),  and  the 
morphological  keys  most  commonly  used  to  identify 
Colombian  Anopheles  (Suarez  et  al.  1988)  predate  the 
description  of  An.  calderoni.  Following  the  detection  of 
atypical  An.  punctimacula  s.l.  specimens  from  various 
localities  in  the  valley  of  the  Rio  Cauca,  progeny  broods 
were  obtained  to  facilitate  a  thoroughly  integrated  sys¬ 
tematic  assessment  of  these  Colombian  variants.  Identi¬ 
fications  of  archive  specimens  of  An.  punctimacula  s.l. 
in  the  Entomological  Museum,  Universidad  del  Valle, 
Colombia  (MUSENUV)  were  revised,  alongside  wild- 
caught  material  from  Ecuador. 
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MATERIALS  AND  METHODS 

Specimens  -  Wild-caught  adult  females  of  An.  punc- 
timacula  sd.  were  collected  in  August  2007,  between  18 
pm-21  pm,  either  resting  or  in  human  landing  catches, 
at  Corporacion  Autonoma  Regional  del  Valle  del  Cauca 
in  Laguna  de  Sonso,  Buga,  Valle  del  Cauca,  Colombia 
(03°52’32”N,  76°20’53”W)  at  an  elevation  of  960  m.  Fe¬ 
males  were  induced  to  oviposit  according  to  the  standard 
protocols  of  Belkin  et  al.  (1965)  and  once  eggs  were  laid, 
females  were  pinned  and  stored  in  the  MUSENUV  collec¬ 
tions  and  labelled  with  the  same  code  as  the  correspond¬ 
ing  progeny.  Larvae  were  fed  on  a  mix  of  brewer’s  yeast 
and  rabbit  food.  Some  4th  instar  larvae  were  separated, 
killed  by  boiling  water  and  stored  in  acetic  acid,  glyc¬ 
erol  and  alcohol  fixative  (10  parts  60%  ethyl  alcohol  with 
1  part  glacial  acetic  acid  and  1  part  glycerol)  for  future 
morphological  examination.  Larval  and  pupal  exuviae  of 
link-reared  specimens  were  stored  in  75%  ethanol  prior  to 
slide-mounting  following  the  protocols  in  Gonzalez  and 
Carrejo  (2007).  Emergent  adults  were  allowed  to  feed  on 
10%  sucrose  solution  and  were  later  killed  by  freezing. 
Adults  were  pinned  and  labelled  with  the  same  code  as 
the  mother  for  further  integrated  studies. 

Morphological  assessment  -  Morphological  identi¬ 
fications  were  carried  out  using  the  diagnostic  charac¬ 
ters  for  An.  punctimacula  s.l.  (Wilkerson  1990)  and  the 
original  species  description  of  An.  calderoni  (Wilkerson 
1991).  Wing  character  nomenclature  used  herein  follows 
that  of  Wilkerson  and  Peyton  (1990)  (Fig.  1,  Tables  I-III). 
Measurements  were  taken  using  a  micrometer  mounted 
in  a  Nikon  SMZ-645  Stereo  Zoom  Microscope.  For  each 
female  (n  =  79),  variation  in  the  overlap  of  the  postsub- 
costal  dark  spot  with  the  preapical  dark  spot  (Table  I), 
the  respective  lengths  of  the  costal  wing  spots  (Table  II) 
and  the  ratio  of  the  length  of  the  costal  wing  spot  rela¬ 
tive  to  the  total  wing  length  (measured  from  the  humeral 
cross-vein)  (Table  III)  were  calculated.  The  presence  of 
key  diagnostic  characters  of  An.  calderoni  larvae  -  num¬ 
bers  of  branches  on  seta  3-C  and  the  positioning  of  seta 
1-X  in  relation  to  the  saddle  (Wilkerson  1991)  -  were  as¬ 
sessed  in  all  mounted  larvae  (n=  144)  (Table  IV).  Verifi¬ 
cation  of  specimens  held  in  the  archive  collections  of  the 
MUSENUV  and  which  had  been  identified  as  An.  punc¬ 
timacula  were  also  reassessed  (Supplementary  data). 

Molecular  analysis  -  Full  locality  details  of  speci¬ 
mens  used  in  the  molecular  part  of  this  study  are  given 
in  Table  V.  DNA  extraction,  polymerase  chain  reaction 
(PCR)  amplification  and  product  clean-up  took  place  on 
board  The  Scholarship  (M/V  Oceanic  II)  as  part  of  the 
activities  of  the  Mosquito  Barcoding  Initiative  (MBI). 
DNA  was  extracted  from  the  abdomens  of  pinned  speci¬ 
mens  using  the  QIAGEN®  BioSprint  robotic  DNA  ex¬ 
traction  machine  and  the  BioSprint  96  DNA  Tissue  Kit 
(QIAGEN®,  Crawley,  UK),  following  the  manufacturers 
instructions.  The  universal  LCO  and  HCO  barcoding 
primers  of  Folmer  et  al.  (1994)  were  used  to  amplify 
the  cytochrome  c  oxidase  subunit  I  gene  (COI)  region 
(658  bp,  less  primers).  The  PCR  comprised  1  |iL  template 
DNA,  1  pL  10  x  NH4  buffer,  0.5  pL  dNTPs  at  2.5  mM, 


0.3  pL  each  primer  at  10  pM,  0.4  pL  MgCL,  at  50  mM 
and  0.2  pL  of  Taq  polymerase  (BioLine,  London,  Eng¬ 
land),  made  up  to  10  pL  with  ddH,0.  The  PCR  reaction 
comprised  initial  denaturing  at  95°C  for  5  min,  then 
34  cycles  of  95°C  for  30  s,  48°C  for  30  min  and  72°C  for 
45  s,  followed  by  a  5-min  extension  at  72°C  and  a  10°C 
hold.  Two  microliter  of  each  PCR  product  was  visual¬ 
ized  on  1%  agarose  gels,  containing  0.5  mg/mL  of  ethid- 
ium  bromide.  Amplification  of  the  nuclear  ITS2  rDNA 
region  was  achieved  using  the  5.8SF  and  28SR  primers 
of  Collins  and  Paskewitz  (1996),  following  the  protocol 
listed  in  Linton  et  al.  (2001). 

PCR  products  were  purified  using  ExoSAP-IT®  (USB, 
Cleveland,  Ohio),  following  the  manufacturers  instruc¬ 
tions.  Bi-directional  DNA  sequences  were  generated  ei¬ 
ther  in  the  Sequencing  Facility  of  the  Natural  History  Mu¬ 
seum,  London,  or  at  the  Museum  Support  Center  of  the 
Smithsonian  Institution,  Maryland,  using  the  Big  Dye®’ 
Terminator  Kit  (PE  Applied  Biosystems,  Warrington, 
England)  and  run  on  ABI  3730  automated  sequencers 
(PE  Applied  Biosystems).  Sequences  were  edited  using 


Subcostal  area 
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Fig.  1:  schematic  diagram  of  the  wing  of  Anopheles  calderoni.  A: 
the  position  of  postsubcostal  pale  spot  (POSCP)  with  relation  to  its 
position  on  the  first  radial  vein  (Rj);  B:  variation  in  the  position  of 
the  POSCP  with  regard  to  its  reportedly  standard  position  on  R  . 
PD:  preapical  dark;  POSCD:  postsubcostal  dark.  PRSCD:  presubcos- 
tal  dark;  PRSCP:  presubcostal  pale;  SCD:  subcostal  dark. 
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TABLE  I 

Variation  in  the  extent  of  overlap  (in  mm)  of  the  distal  part  of  postsubcostal  dark  spot  with  the  proximal  part  of  the  preapical 
dark  spot  on  vein  first  radial  vein  from  79  female  progeny  of  iso-female  lines  of  Anopheles  calderoni 
generated  from  females  collected  in  Valle  de  Cauca,  Colombia 


Iso  family 

Iso  family  number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1° 

0.020 

-0.010 

0.080 

0.010 

0.060 

0.040 

0.050 

-0.040 

0.100 

0.020 

2 

0.040 

-0.040 

0.110 

0.026 

0.090 

0.080 

-0.060 

0.000 

0.040 

0.060 

3 

0.030 

0.000 

0.050 

0.054 

0.060 

0.000 

0.000 

0.010 

0.036 

0.004 

4“ 

0.026 

-0.026 

0.074 

0.020 

0.060 

0.070 

0.070 

0.000 

0.050 

0.000 

5° 

0.000 

-0.020 

0.040 

0.040 

0.090 

0.020 

0.110 

0.000 

0.006 

0.070 

6 

- 

-0.020 

0.070 

0.060 

0.054 

0.030 

0.000 

-0.040 

0.026 

0.050 

7 

- 

-0.034 

0.080 

- 

0.130 

0.000 

0.020 

-0.030 

-0.010 

0.020 

8 

- 

-0.040 

0.050 

- 

- 

0.020 

0.000 

- 

- 

0.020 

9 

- 

-0.030 

0.026 

- 

- 

- 

-0.040 

- 

- 

0.000 

10 

- 

-0.050 

0.050 

- 

- 

- 

-0.030 

- 

- 

- 

Mean 

0.023 

-0.024 

0.063 

0.035 

0.078 

0.033 

0.012 

-0.014 

0.035 

0.027 

SD 

±0.015 

±0.015 

±  0.024 

±  0.020 

±  0.028 

±  0.030 

±  0.052 

±0.021 

±  0.035 

±  0.026 

a:  DNA-typed  isofamily;  SD:  standard  deviation.  Negative  or  zero  values  signify  lack  of  overlap. 


Sequencher™  version  4.8  (Genes  Codes  Corporation, 
Ann  Arbor,  Michigan)  and  aligned  in  CLUSTALX  (Jean- 
mougin  et  al.  1998)  (Fig.  2).  Similarity  of  the  sequences 
generated  in  this  study  was  compared  with  those  avail¬ 
able  in  GenBank using  Blast  (blast.ncbi.nlm.nih.gov/)  and 
sequences  were  analyzed  according  to  standard  barcod¬ 
ing  methodology  (Hebert  et  al.  2003).  Sequence  statistics, 
calculation  of  pairwise  distance  parameters  using  Kimura 
2-parameter  algorithm  (Kimura  1980)  and  generation  of 
the  bootstrapped  neighbor-joining  tree  (10,000  replicates) 
were  carried  out  in  MEGA  v.4.0  (Tamura  et  al.  2007). 

Details  of  isofamily  specimens  of  An.  calderoni  se¬ 
quenced  from  Colombia  (n  =  7)  and  wild-caught  speci¬ 
mens  from  Ecuador  (n  =  5)  in  this  study  are  given  in 
Table  V.  In  addition,  one  specimen  each  of  the  closely 
related  taxa  An.  punctimacula  [Gamboa,  Panama 
(09°06’N,  79°41’W);  colls.:  J  Loaiza,  Y  Linton,  F  Ruiz 
and  students  of  The  Scholarship  LLC]  and  An.  malefac¬ 
tor  [Yaviza,  Panama  (08°11’N,  77°41’W);  coll.:  J  Loaiza] 
were  included  in  the  analysis,  but  no  specimens  of  An. 
guarao  were  available  for  inclusion.  The  COI  electro- 
pherograms  and  specimen  details  (including  exact  lo¬ 
calities  details  with  georeferences,  details  of  voucher 
type  and  housing,  GenBank  accessions)  are  available  in 
the  public  records  of  the  MBI  section  of  the  Barcode  of 
Life  Data  Systems  v.  2.5  (boldsystems.org)  (Ratnasing- 
ham  &  Hebert  2007)  and  appear  in  GenBank  as  barcode 
red  flag  data,  indicating  their  high  quality  and  voucher 
standards  under  the  following  accessions:  An.  calderoni 
(Colombia  HQ642964-70,  Ecuador  HQ642971-74),  An. 
punctimacula  (Panama  HQ622626)  and  An.  malefactor 
(Panama  HQ622625).  The  ITS2  sequences  of  An.  calde¬ 
roni  generated  in  this  study  are  available  in  GenBank 
under  accessions  numbers  HQ622618-HQ622622. 


RESULTS 

Morphological  assessment  of  Colombian  An.  calderoni 

Adults  -  Nearly  all  the  female  progeny  from  families 
reared  for  this  study  exhibited  characters  concordant  with 
those  given  in  the  original  description  of  An.  calderoni 
as  follow:  upper  mesanepimeron  with  pale  scales,  wings 
with  pale  yellow  scales  mixed  with  white  scales,  distal 
postsubcostal  pale  spot  on  the  costa  and  the  postsubcostal 
pale  spot  on  the  first  radial  (R^  not  contiguous,  frequently 
separated  by  dark  scales  at  their  respective  ends,  i.e.  over¬ 
lap  of  scales  of  the  distal  part  of  the  postsubcostal  dark 
spot  with  the  proximal  part  of  the  preapical  dark  on  vein 
R,  (Fig.  1A,  B);  scales  on  vein  R[  between  the  proximal 
part  of  the  presubcostal  pale  and  the  distal  postsubcostal 
pale  generally  dark  (Fig.  IB),  i.e.  R  in  the  subcostal  area 
is  usually  all  dark  except  for  pale  spots  at  the  ends. 

Measures  of  the  overlap  of  scales  of  the  distal  part  of 
the  postsubcostal  dark  spot  with  the  proximal  part  of  the 
preapical  dark  spot  of  R[  were  found  to  be  heterogeneous 


SSsSSSSissSIs 
Hap1otypel(COL=6)  tgcggcaaaacta 

Hap1otype2(coL=l)  ....a . 

Haplotype3(EC=2)  c. . .at. .g. .eg 
Haplotype4(EC=l)  c. . .atg.g. .eg 
Hapl otype5 (EC=1)  ca . aa . gg . . ttg 
Hapl otype6(EC=l5  cataa . g . . getg 

Fig.  2:  alignment  of  unique  cytochrome  c  oxidase  subunit  I  gene  hap- 
lotypes  of  13  Anopheles  calderoni  from  Colombia  (COL)  (n  =  7)  and 
Ecuador  (EC)  (n  =  6).  Numbers  above  bases  indicate  position  of  vari¬ 
able  bases  along  the  658  bp  fragment  (709  bp  with  primers).  Points 
indicate  bases  identical  to  haplotype  1. 
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TABLE  III 


Measurements  of  the  costal  wing  spots  normalized  for  wing  length  (humeral  cross  vein  to  apex)  from  79 
females  from  10  isofamilies  of  Anopheles  calderoni  generated  from  females  collected  in  Valle  de  Cauca, 
Colombia,  compared  to  those  given  in  the  original  description  of  An.  calderoni  (Wilkerson  1991) 


Wing  spot 

Iso  families  (Colombia) 

Wilkerson  (1991) 

Mean 

SD  (±) 

Range 

Mean 

SD  (±) 

Range 

Humeral  pale 

0.01 

0.010 

0.01-0.03 

0.02 

0.010 

0.01-0.03 

Humeral  dark 

0.07 

0.010 

0.04-0.10 

0.09 

0.020 

0.03-0.12 

Presector  pale 

0.01 

0.000 

0.01-0.02 

0.02 

0.010 

0.01-0.03 

Presector  dark 

0.11 

0.020 

0.08-0.15 

0.09 

0.020 

0.05-0.14 

Sector  pale 

0.07 

0.020 

0.03-0.10 

0.10 

0.020 

0.06-0.14 

Sector  pale  accessory 

0.04 

0.010 

0.02-0.07 

0.06 

0.020 

0.04-0.10 

Sector  dark 

0.18 

0.020 

0.02-0.23 

0.15 

0.020 

0.11-0.18 

Area  subcostal 

0.24 

0.020 

0.20-0.29 

0.28 

0.020 

0.24-0.31 

Presubcostal  dark 

0.04 

0.010 

0.01-0.06 

0.04 

0.010 

0.02-0.06 

Presubcostal  pale  (proximal) 

0.02 

0.010 

0.01-0.05 

0.02 

0.010 

0.01-0.06 

Presubcostal  pale  (distal) 

0.02 

0.000 

0.01-0.04 

0.02 

0.010 

0.01-0.03 

Postsubcostal  dark 

0.10 

0.020 

0.06-0.14 

0.10 

0.010 

0.07-0.12 

Postsubcostal  pale  (proximal) 

0.02 

0.000 

0.01-0.03 

0.02 

0.010 

0.01-0.04 

Postsubcostal  pale  (distal) 

0.02 

0.010 

0.02-0.04 

0.04 

0.010 

0.02-0.06 

Subcostal  dark 

0.03 

0.010 

0.01-0.05 

0.02 

0.001 

0.01-0.04 

Preapical  dark 

0.15 

0.020 

0.12-0.20 

0.13 

0.016 

0.10-0.18 

Preapical  pale 

0.10 

0.010 

0.07-0.12 

0.11 

0.010 

0.09-0.12 

Preapical  dark  accessory 

0.05 

0.010 

0.01-0.07 

0.05 

0.010 

0.04-0.07 

Apical  dark 

0.03 

0.010 

0.01-0.04 

0.02 

0.010 

0.01-0.03 

SD:  standard  deviation. 


TABLE  IV 

Variation  in  numbers  of  branches  of  seta  3-C  and  the  positioning  of  seta  1-X  with  respect  to  the  saddle  in  144  specimens  from 
10  isofamilies  of  Anopheles  calderoni  generated  from  females  collected  in  Valle  de  Cauca,  Colombia 


Isofamily 

Individuals 

n 

Number  of  branches  on 

seta  3-C 

Position  of 

seta  1-X 

Mean 

Min 

Max 

Mode 

n 

Saddle 

% 

Border 

% 

1" 

6 

5.5 

4 

7 

4-7 

3 

100 

- 

2 

6 

6.2 

4 

7 

4-8 

3 

100 

- 

3 

26 

5.9 

3 

9 

4 

13 

69.2 

23.1 

4“ 

26 

4.5 

3 

7 

4 

12 

91.7 

8.3 

5“ 

13 

4.5 

2 

6 

4-6 

6 

100 

- 

6 

5 

6.6 

5 

7 

5 

3 

100 

- 

7 

19 

5.4 

2 

7 

6 

10 

70 

30 

8 

12 

5.8 

4 

8 

6 

5 

60 

40 

10 

15 

4.7 

1 

6 

6 

7 

85.7 

14.3 

11 

5 

7.4 

1 

11 

- 

2 

50 

50 

12 

4 

5.5 

5 

6 

5-6 

2 

100 

- 

Total 

137 

5.4 

1 

11 

6 

66 

81.8 

16.7 

a:  DNA-typed  isofamilies. 
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among  families  (p  <  0.0001)  and  variable  even  among  fe¬ 
male  siblings  (Table  I).  This  value  is  positive  when  there  is 
overlap  in  the  spots,  but  can  also  have  a  value  of  zero  when 
the  ends  of  the  spot  are  immediately  adjacent  and  a  nega¬ 
tive  value  when  the  spots  are  clearly  separated  (Fig.  IB). 
Overall,  eight  of  the  10  isofamilies  examined  exhibited 
overlapping  spots  (0.003-0.024  mm),  however  many  of  the 
individuals  in  isofamily  2  and  some  from  isofamilies  7-9 
showed  either  no  overlap  but  immediately  adjacent  spots  or 
clear  separation  of  the  spots  (Table  I).  An  even  alignment 
of  the  postsubcostal  spot  on  C  and  Rp  reportedly  typical 
of  all  other  Arribalzagia  Series  species  (Suarez  et  al.  1988, 
Wilkerson  1990, 1991,  Gonzalez  &  Carrejo  2007),  was  not 
observed  in  the  An.  calderoni  specimens  examined  here. 
Considering  the  above  observations,  it  is  necessary  to  add 
a  clarification  of  the  original  description  of  this  character 
as  follows:  “distal  postsubcostal  pale  spot  on  vein  C  gener¬ 
ally  not  contiguous  with  the  postsubcostal  pale  on  vein  Rp 
separated  (generally)  or  not  by  dark  scales  of  the  distal  part 
of  the  postsubcostal  dark  and  the  proximal  part  of  preapi- 
cal  dark  on  R  For  those  specimens  that  agree  with  the  re¬ 
maining  characters  yet  where  the  aforementioned  overlap 
does  not  occur,  it  is  possible  complement  species  identifi¬ 
cation  by  examination  of  the  scales  between  the  proximal 
presubcostal  pale  and  the  distal  postsubcostal  pale  spots 
located  on  Rr  These  are  usually  all  dark  in  An.  calderoni, 
whereas  they  are  predominantly  pale  in  An.  punctimacula. 
An.  malefactor  and  An.  guarao  (Wilkerson  1991). 

Analysis  of  costal  spots  within  the  isofamilies  re¬ 
vealed  that  some  specimens  did  not  comply  with  a  nor¬ 
mal  distribution,  even  after  data  transformation  using 
wing  length;  hence  a  non-parametric  Kruskal-Wallis 
test  was  used  to  compare  the  averages  of  these  measure¬ 
ments  among  isofamilies. 

Statistically  significant  differences  were  found  when 
comparing  lengths  of  the  subcostal  area  (Fig.  3)  as  a 
whole  among  the  isofamilies;  yet,  the  relative  size  com¬ 
parisons  of  individual  spots  within  this  area,  i.e.  the  size 
of  the  most  distal  presubcostal  pale  spot,  and  the  proxi¬ 
mal  and  distal  postsubcostal  pale  spots,  did  not  statisti¬ 
cally  differ  amongst  isofamilies  (Table  II).  Furthermore, 
the  average  values  among  isofamilies  (Table  III)  were 
very  similar  to  those  presented  in  the  original  descrip¬ 
tion  (Wilkerson  1991). 

The  band  or  spot  on  hindtarsomere  5  was  generally 
found  not  to  be  as  pale  as  in  the  original  description 
(Wilkerson  1991).  Instead  it  was  observed  as  a  darker 
band,  sometimes  grayish  or  pale  brown  and  was  not  well 
defined  in  91.1%  of  the  progeny  specimens  examined. 
This  darker  band  was  also  noted  in  94.6%  of  specimens 
collected  in  human  landing  catches  from  the  same  local¬ 
ity  (Valle  de  Cauca,  n  =  37). 

The  original  description  of  An.  calderoni  makes  no 
reference  to  variation  in  the  pale  and  dark  spots  of  the 
4th  hindtarsomere  4,  which  are  otherwise  widely  em¬ 
ployed  as  diagnostic  characters  in  several  morphologi¬ 
cal  keys  for  species  of  the  Arribalzagia  Series  (Gorham 
et  al.  1973,  Wilkerson  et  al.  1990,  Gonzalez  &  Carrejo 
2007).  There  are  generally  five  spots  of  varying  length 
on  hindtarsomere  4:  two  dark  and  three  pale  (basal,  api¬ 
cal  and  central).  Spots  reportedly  differ  in  length,  with 


the  two  pale  spots  (basal,  apical)  being  smaller  than  the 
dark  spots.  In  91.6%  of  Colombian  An.  calderoni  exam¬ 
ined  (n  =  79),  the  central  pale  spot  separated  the  two  dark 
spots;  whereas  in  the  remaining  specimens  the  central 
pale  spot  was  not  clear,  allowing  the  two  dark  spots  to 
merge  and  appear  more  like  a  band.  The  central  pale  spot 
was  found  to  be  of  variable  width,  and  was  mostly  light 
grey  in  color.  This  character  was  even  found  to  be  vari¬ 
able  between  legs  of  the  same  individual  therefore  we 
discourage  the  further  use  of  this  character  for  species 
diagnosis  in  the  Arribalzagia  Series. 

Larvae  -  The  number  of  branches  in  seta  3-C  and  the 
position  of  seta  1-X  were  assessed  in  all  144  specimens 
from  10  isofamilies  (Table  IV).  That  seta  1-X  is  not  pres¬ 
ent  on  the  saddle  is  a  shared  characteristic  of  both  An. 
punctimacula  and  An.  calderoni.  In  98.2%  of  progeny  ex¬ 
amined,  seta  1-X  was  not  found  on  the  saddle  itself,  with 
81.5%  attached  outside  this  structure  altogether  and  16.7% 
placed  on  its  border  (Table  IV).  Most  commonly,  An.  cal¬ 
deroni  larvae  exhibited  4-7  branches  on  seta  3-C  (range 
1-11;  mean  5.4;  mode  6).  The  mean  number  of  branches 
is  slightly  higher  than  the  4.2  reported  in  the  original  de¬ 
scription  (Wilkerson  1991),  most  probably  due  to  increased 
sample  size.  An.  punctimacula  differs  significantly  in  hav¬ 
ing  an  average  of  11.5  branches  hence  the  combination  of 
these  two  characters  allows  morphological  differentiation 
of  larval  An.  calderoni  and  An.  punctimacula. 


80"0'0"W  70’0'0-W  60*0'0"W 


Fig.  3:  map  indicating  the  known  distribution  of  Anopheles  calderoni 
in  Colombia,  Ecuador,  Peru  and  Venezuela,  based  on  morphological 
and  molecular  verification,  following  this  study  and  those  of  Wilker¬ 
son  (1991)  and  Rubio-Palis  and  Moreno  (2003). 
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DISCUSSION 

Revision  of  MUSENUV ^ archive  specimens  -  Retrospec¬ 
tive  identification  of  archive  specimens  of  An.  punctimac- 
nla  and  An.  malefactor  held  in  the  MUSENUV  revealed  a 
significant  proportion  of  specimens  were  infact  An.  calde- 
roni.  Specimens  of  An.  calderoni  were  detected  in  collec¬ 
tions  from  1948-2009,  in  35  localities  in  20  municipalities, 
representing  12  departments  of  Colombia  (Supplemen¬ 
tary  data).  Although  the  distribution  of  An.  calderoni  in 
Colombia  appears  predominantly  coastal  (both  on  the  Pa¬ 
cific  and  Atlantic  coasts)  (Fig.  3),  some  of  the  departments 
(e.g.  Tolima)  are  far  inland.  Additionally,  several  collec¬ 
tion  sites  in  Valle  de  Cauca  and  Quindio  were  over  1000  m 
in  altitude  (Supplementary  data),  considerably  higher  than 
the  250  m  limit  previously  reported  for  early  reports  of  An. 
calderoni  in  Peru  (Wilkerson  1991)  and  200-400  m  report¬ 
ed  in  Venezuela  (Rubio-Palis  &  Moreno  2003).  Rather  our 
data  supports  the  findings  of  Calderon  et  al.  (1995)  who 
reported  An.  calderoni  from  sea  level  to  1080  m  in  Peru, 
suggesting  that  the  distribution  of  An.  calderoni  could  be 
considerably  more  extensive  than  at  first  thought. 

Molecular  characterization  of  An.  calderoni  - 
ITS2  rDNA  sequences  generated  from  two  specimens 
[PE400(2)-29,  PE400(2)-31]  belonging  to  the  type  series 
of  An.  calderoni  from  Salitral  (Wilkerson  1991)  were 
identical  to  those  from  three  specimens  from  Ecuador 
(Table  V),  thus  verifying  the  new  country  record  of  An. 
calderoni  in  Ecuador.  Correlation  of  the  658  bp  COI 
DNA  barcode  of  one  of  the  verified  Ecuadorian  An.  cal¬ 
deroni  specimens  (EC229-1)  allowed  the  confirmation 
the  conspecific  identity  of  additional  wild-caught  adults 
from  two  other  sites  in  Ecuador  (Guayas,  n  =  4,  Pichin- 
cha,  n  =  1)  and  progeny  broods  collected  in  Valle  de  Cau¬ 
ca,  Colombia  (n  =  7  from  isofamilies  1,  4  and  5)  (Table 
V),  used  in  the  morphological  assessment  (Table  V). 

In  total,  DNA  barcodes  were  generated  for  12  speci¬ 
mens  (Colombia,  n  =  7,  Ecuador,  n  =  5),  which  varied  in 
AT  richness  from  66-66.5%.  Six  unique  haplotypes  were 
determined,  as  follows:  haplotypes  1  (n  =  6)  and  2  (n  =  1) 
(both  Valle  de  Cauca,  Colombia),  haplotype  3  [Guayas  (n 
=  1)  and  Pichincha  (n  =  1),  Ecuador]  and  haplotypes  4-6  (n 
=  1,  each  one  from  Guayas,  Ecuador)  (Figs  2,  3).  Within 
Colombia,  haplotypes  1  and  2  varied  by  only  one  base 
(A/G  at  base  319)  (Fig.  2).  Overall,  13  (1.97%)  variable 
bases  were  noted,  10  of  which  were  parsimony  informa¬ 
tive.  Mean  pairwise  (K2P)  sequence  diversity  within  An. 
calderoni  was  0.70%,  well  within  values  of  intraspecific 
variation  reported  for  the  barcoding  region  in  other  mos¬ 
quito  species  (Cywinska  et  al.  2006,  Kumar  et  al.  2007, 
Azari-Hamidian  et  al.  2010).  Flaplotypes  were  unique 
by  country,  with  mean  intra-country  pairwise  distances 
of  0.04%  and  0.77%  within  Colombian  and  Ecuadorian 
specimens,  respectively,  and  an  overall  mean  inter-coun¬ 
try  sequence  divergence  of  1.03%.  Amino  acid  translation 
of  the  COI  sequences  revealed  that  all  nucleotide  changes 
were  synonymous,  suggesting  that  these  represent  popu¬ 
lation  level  differences  only. 

Differentiation  of  An.  calderoni,  An.  punctimacula, 
An.  malefactor  and  An.  guarao  -  Given  that  An.  calde¬ 
roni  has  often  been  misidentified  as  An.  punctimacula  in 


Colombia  (herein),  we  set  out  to  assess  the  utility  of  DNA 
barcodes  for  species  differentiation  between  these  species 
and  other  closely  related  taxa  -  An.  malefactor  and  An. 
guarao.  No  specimens  of  An.  guarao  were  available  for 
molecular  comparison,  but  the  morphological  distinctive¬ 
ness  of  this  species  makes  it  unlikely  to  be  confused  with 
An.  calderoni.  The  robust  boot  strap  values  illustrated  on 
the  neighbor-joining  tree  generated  using  the  COI  bar¬ 
code  data  of  An.  calderoni,  An.  punctimacula  and  An. 
malefactor  confirmed  the  separate  species  status  of  these 
taxa  (Fig.  4).  In  comparison  to  the  0.70%  interspecific 
variation  noted  within  An.  calderoni,  interspecific  com¬ 
parisons  with  those  of  Panamanian  An.  punctimacula  and 
An.  malefactor  generated  in  this  study,  revealed  5.8%  and 
10.9%  variation,  respectively  (Fig.  4).  COI  barcodes  of 
An.  malefactor  (67.7%  AT)  were  most  similar  to  those  of 
An.  punctimacula  (67.1%  AT),  differing  by  9.6%. 

Neither  An.  calderoni  nor  An.  malefactor  axe  current¬ 
ly  represented  in  GenBank  and  thus  these  sequences  also 
comprise  the  first  DNA  sequences  reported  for  these  spe¬ 
cies.  Although  the  GenBank  entry  AF417719  of  Sallum 
et  al.  (2002)  partially  overlapped  our  An.  punctimacula 
sequence  by  93  bases,  the  sequences  generated  herein 
represent  the  first  published  DNA  barcodes  for  that  spe¬ 
cies.  This  study  shows  that  DNA  barcodes  will  aid  cor¬ 
rect  identification  of  An.  calderoni,  An.  malefactor  and 
An.  punctimacula  in  future  studies  and  will  allow  the 
vector  status  of  An.  calderoni  and  An.  malefactor  to  be 
more  accurately  assessed. 

DNA  verification  of  Ecuadorian  and  Colombian 
specimens  as  An.  calderoni  along  with  the  morphologi¬ 
cal  reassessment  of  archived  MUSENUV  specimens 
significantly  adds  to  the  knowledge  of  the  distribution  of 
An.  calderoni  in  Latin  America  (Fig.  3).  Given  that  the 
known  distribution  of  An.  calderoni  now  includes  Ecua¬ 
dor  and  Colombia,  as  well  as  Peru  and  Venezuela,  and 
its  prior  taxonomic  confusion  with  An.  punctimacula 
and  An.  malefactor,  it  is  now  prudent  to  carry  out  further 
collections  to  determine  the  true  distributions  and  spe¬ 
cific  roles  in  malaria  transmission  of  these  three  species 
across  the  Neotropical  Region. 


Fig.  4:  neighbor-joining  tree  of  Anopheles  calderoni ,  Anopheles  punc¬ 
timacula  and  Anopheles  malefactor  based  on  658  bp  of  the  mtDNA 
cytochrome  c  oxidase  subunit  I  gene  gene.  The  tree  was  constructed 
using  10,000  replicates  of  pairwise  genetic  distances  generated  using 
the  Kimura  2-parameter  distance  algorithm.  Bootstrap  confidence 
limits  are  shown  on  tree  labels. 
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':  specimens  from  the  type  series  of  An.  calderoni;  MUSENUV:  Entomological  Museum,  Universidad  del  Valle,  Colombia. 
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Collection  details  of  specimens  of  Anopheles  calderoni  retrospectively  identified  in  the  collections 
of  the  Entomology  Museum  of  the  Universidad  del  Valle,  Colombia 


Elevation 

Collection 

Department 

Municipality 

Locality 

(m) 

date 

Collector 

Gomez  Plata 

El  Brasil 

244 

10  Mayl984 

SEM 

Antioquia 

Puerto  Nare 

El  Pescado 

270 

4  Oct  1983 

SEM 

Santafe  de  Antioquia 

Paso  Real 

666 

17  Mar  1984 

SEM 

Bolivar 

Achi 

Platanal 

23 

26  Jun  1984 

SEM 

Caldas 

La  Dorada 

Puerto  Tejada 

Guarinocito 

Perico  Negro 

200 

990 

unknown 

unknown 

SEM 

SEM 

Cauca 

Villarica 

Agua  Azul 

990 

18  Feb  2005 

1  Mar  2006 

Curso  de  E.  Medica 

Choco 

Riosucio 

Cano  Pavas 

4 

5  Jun  1981 

SEM 

Guajira 

Riohacha 

Las  Flores 

9 

22  Jun  1986 

H  Ramirez  (SEM) 

Magdalena 

Santa  Marta 

Guachaca 

18 

2  May  1983 

Derrueco  (SEM) 

Brazo  Rio  Mira 

15 

7  Apr  1986 

SEM 

Bucheli 

10 

9  Jun  1987 

Cortes  (SEM) 

Chajal 

9 

7  Jun  1987 

Cortes  (SEM) 

17  Jun  1986 

Chilvi-Linea 

15 

22  Jul  1986 

21  Jul  1986 

Cortes  (SEM) 

Narino 

23  Jul  1986 

El  Descalzadero 

15 

28  Feb  1984 

SEM 

23  Jun  1987 

Ing.  del  Carmen 

15 

3  Jul  1987 

10  Jul  1987 

16  Jul  1987 

Cortes  (SEM) 

Ing.  del  Guadual 

15 

18  Jun  1987 

Cortes  (SEM) 

Nueva  Union  2 

15 

1  Apr  1986 

Castillo  (SEM) 

Norte  de 

Santander 

Tibu 

Puerto  Barco 

Isla  La  Ceiba 

56 

35 

17  Jun  1986 

26  Jun  1987 

24  Sep  1986 

Pabon  (SEM) 

SEM 

Cano  Castillo 

60 

Sep  1983 

SEM 

Quindio 

Quimbaya 

La  Montana  del  Ocaso 

1100 

13  Nov  2003 

P  Lopez  &  I  Valderrama 

To  lima 

Ataco 

Bocas  de  S.  Pedro 

613 

30  Nov  1981 

D  Tique  (SEM) 
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Elevation 

Collection 

Department  Municipality 

Locality 

(m) 

date 

Collector 

R  Gonzalez,  N  Carrejo, 

7  Jun  2007 

DM  Lucumi,  CX 

Laguna  de  Sonso 

960 

9  Jul  2007 

Restrepo, 

Buga 

May  2009 

DM  Montoya  &  C 

Salas 

El  Vinculo 

970 

2008 

LM  Martinez 

Chambimbal,  El  Cedral 

928 

May  2009 

D  M  Montoya 

Bugalagrande 

El  Pital  (Humedal) 

926 

29  Apr  2009 

DM  Montoya 

Candelaria 

El  Cabuyal 

El  Otono 

1005 

1005 

May  2006 

7  May  2007 

May  2008 

MH  Cuellar  &  LM 

Martinez 

DM  Cundumi  &  CX 

Restrepo 

Valle  del  Cauca 

Chocosito 

1045 

15  May 

2007 

R  Gonzalez 

Florida 

May  1948 

Florida  (Cab) 

1030 

20  Aug 

A  Patino  (SEM) 

1948 

Jamundi 

Urb.  Terranova 

980 

28  Apr  2009 

V  Ceron 

Ciat 

962 

Oct  2007 

LM  Martinez 

15  May 

Palmira 

San  Isidro-El  Bolo 

980 

2005 

24  Nov 

D  Gomez  &  LM 

Martinez 

2007 

Tulua 

Mate  Guadua 

1113 

28  Apr  2006 

H  Cadena,  NJ  Mina  & 

C  Robledo 

SEM:  Servicio  de  Erradicacion  de  la  Malaria,  Colombia. 


